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ABSTRACT 

We present a preliminary measurement of the rate of gluon splitting into bottom quarks, 
g — > bb, in hadronic Z° decays collected by SLD between 1996 and 1998. The analysis 
was performed by looking for secondary bottom production in 4-jet events of any 
primary flavor. 4-jet events were identified, and a topological vertex-mass technique 
was applied to each jet in order to identify b or b jets. The upgraded CCD based 
vertex detector gives very high 5-tagging efficiency, especially for B hadrons of the 
low energies typical of this process. The two most nearly collinear bjb jets were tagged 
as originating from g — ► bb. We measured the rate of secondary bjb production per 
hadronic event, g b i, to be (3.07 ± 0.71(stat.) ± 0.66(syst.)) x 10 -3 (preliminary). 
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1. Introduction 



The process of the splitting of a gluon into a heavy-quark pair is one of the elemen- 
tary processes in QCD but is poorly known, both theoretically and experimentally. 

The rate g bb is defined as the fraction of hadronic events in which a gluon splits 
into a bb pair, e + e~ — > qqg — > qqbb. The value of g bb is an infrared finite quantity, 
because the 6-quark mass provides a natural cutoff, hence it can be safely computed in 
the framework of perturbative QCD However the rate is sensitive to as and to the 
6-quark mass, which results in a substantial theoretical uncertainty in the calculation of 
g bb . The limited accuracy of the g bb prediction is one of the main sources of uncertainty 
in the measurement of the partial decay width R b = T(Z° — > bb)/T(Z° — > qq) |fj, [|. 
In addition, about 50% of the B hadrons produced at the Tevatron are due to the 
gluon splitting process, and a larger fraction is expected to contribute at the LHC. 
A better knowledge of this process can improve theoretical predictions of heavy-flavor 
production at such colliders. 

This measurement is difficult experimentally. The cross section for g — > bb is very 
small even at Z° energies, since the gluon must have sufficient mass to produce the 
bottom-quark pair. There are huge backgrounds from Z° — > bb whose magnitude is 
about a hundred times larger than the Z° — > qqg — > qqbb process. Moreover the B 
hadrons from g — ► bb have relatively low energy and short flight distance and are more 
difficult to distinguish using standard vertexing. So far, measurements of g bb have been 
reported by DELPHI and ALEPH g. 

Here we present a new measurement of g bb based on a 400k Z°-decay data sample 
taken in 1996-98 at the Stanford Linear Collider (SLC), with the SLC Large Detector 
(SLD). In this period, Z° decays were collected with an upgraded vertex detector, wider 
acceptance and better impact parameter resolution, thus improving considerably the 
6-tagging performance. 
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2. The SLD Detector 



A description of the SLD is given elsewhere ||. Only the details most relevant to this 
analysis are mentioned here. 

SLD is well-suited for the measurement of (7 — ^ 66 due to two unique features. 
The first is that the SLC, the only linear collider in the world, provides a very small 
and stable beam spot. The SLC interaction point was reconstructed from tracks in 
sets of approximately thirty sequential hadronic Z° decays with an uncertainty of 
only 5/im transverse to the beam axis and 32/zm (for bb events) along the beam axis. 
Second is the upgraded vertex detector (VXD3) 0, a pixel-based CCD vertex detector. 
VXD3 consists of 3 layers with 307M pixels and each layer is only 0.36% of a radiation 
length thick. The measured rip (rz) track impact-parameter resolution approaches 
11/zm (23/im) for high momentum tracks, while multiple scattering contributions are 
40/zm/(j9sin 3//2 9) in both projections (z is the coordinate parallel to the beam axis). 
With these features, topological vertex finding gives excellent 6-tagging efficiency and 
purity. In particular, the efficiency is good even at low B-meson energies, which is 
especially important for detecting g — > bb. 



3. Flavor Tagging 



Topologically reconstructed secondary vertices [0] are used by many analyses at the 
SLD for heavy-quark tagging. To reconstruct the secondary vertices, the space points 
where track density functions overlap are found in 3-dimensions. Only the vertices 
that are significantly displaced from the primary vertex (PV) are considered to be 
possible B- or D-hadron decay vertices. The mass of the secondary vertex is calculated 
using the tracks that are associated with the vertex. Since the heavy-hadron decays 
are frequently accompanied by neutral particles, the reconstructed mass is corrected to 
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account for this fact. By using kinematic information from the vertex flight path and 
the momentum sum of the tracks associated with the secondary vertex, we calculate 
the Py-corrected mass Mp T by adding a minimum amount of missing momentum to 
the invariant mass, as follows: 

VTX + P T 2 + \P T \. 

Here Myrx is the invariant mass of the tracks associated with the reconstructed sec- 
ondary vertex and Pt is the transverse momentum of the charged tracks with respect 
to the B-fTight direction. In this correction, vertexing resolution as well as the PV 
resolution are crucial. Due to the small and stable interaction point at the SLC and 
the excellent vertexing resolution from the SLD CCD Vertex detector, this technique 
has so far only been successfully applied at the SLD. 



4. Monte Carlo and data Samples 

The measurement uses 400k Z° — > hadron events collected between 1996 and 1998 with 
the requirement that the VXD3 was fully operational. 

For the purpose of estimating the efficiency and purity of the g — > bb selection 
procedure, we made use of a detailed Monte-Carlo simulation of the detector. The 
JETSET 7.4 || event generator was used, with parameter values tuned to hadronic 
e + e~ annihilation data [[J, combined with a simulation of B hadron decays tuned to 



T(4S) data |TIJ and a simulation of the SLD based on GEANT 3.21 |TTJ. Inclusive 
distributions of single-particle and event-topology observables in hadronic events were 
found to be well described by the simulations . Uncertainties in the simulation were 
taken into account in the systematic errors (Section |7|). 

Monte-Carlo events are reweighted to take into account current estimates for gluon 
splitting into heavy-quark pairs [|], fI3|| . JETSET with the SLD parameters predicts 
g h l = 0.14% and g c5 = 1.36%, and we reweighted them so that g b i = 0.273% and 



g c5 = 2.58%. A Monte-Carlo production of about 1200k Z — > qq events, 1000k Z — > bb 
events and 480k Z — > cc events are used in order to better evaluate the efficiencies. 

Besides the signal events, hereafter called B, two categories of background events 
exist: 

• Events which do not contain any gluon splitting into heavy flavor at all, hereafter 
called Q events; and 

• Events in which a gluon splits to a charm quark pair, named C events. 



5. Analysis 

The two B hadrons coming from the gluon tend to be produced in a particular topo- 
logical configuration, which allows one to discriminate the signal from background. We 
select g — > bb events as follows: 

• Require 4 jets in the events; 

• Require b tags in two jets selected in a particular configuration; and 

• Apply additional topological selections to improve the signal/background ratio. 
Jets are formed by applying the Durham jet-finding algorithm |14| to energy-flow 



particles with y cut = 0.008, chosen to minimize the statistical error. The overall 4- 
jet rate in the data is (5.976 ± 0.044)%, where the error is statistical only. In the 
Monte-Carlo simulation the 4-jet rate is (5.678 ± 0.002 ± 0.068)% where the first error 
is statistical and the second is due to the uncertainty in the simulation of heavy-quark 
physics. The 4-jet rates for the B, C and Q events predicted by the simulation are 
about 32%, 18% and 5.3%, respectively. The two jets forming the smallest angle in 
the event are considered as candidates for originating from the gluon splitting process 
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g — > bb. The selected jets are labeled as jet 1 and jet 2, where jet 1 is more energetic 
than jet 2. The other two jets in the event are labeled as jets 3 and 4, where jet 3 is 
more energetic than jet 4. 

Jets containing B-hadron decay products are then searched for by making use of 
the information coming from the vertex detector, using the topological vertex method. 
We require both jet 1 and jet 2 to contain a secondary vertex. No tag is applied to jet 
3 and jet 4. After topological vertexing, about 300 events are selected. The selection 
efficiency for g — > bb is expected from the Monte Carlo simulation to be 6.6% while 
the signal/background ratio is 1/5. 67% of the background comes from Z —>■ bb events, 
21% from g — > cc events and remaining 12% from Z — > qq (q ^ b) events. 

In order to improve the signal/background ratio, we use topological information. 
Firstly, many Z° — > bb background events have one 6-jet which was split by the jet- 
finder into 2 jets so that the two found vertices are from different decay products from 
the same B decay. The two vertex axes tend to be collinear. Figure |l] shows the angular 
distribution between vertex axes in jet 1 and jet 2. Half of the bb background peaks at 
cos 9 12 ~ 1. In order to remove bb events, we require —0.2 < cos6 , 12 < 0.96. 

Secondly, the variable | cos 01234!, where 01234 is the angle between the plane n 12 
formed by jets 1 and 2 and the plane II34 by jets 3 and 4, is used to suppress the bb 
background. Figure |^ shows the distribution of | cos 01234 1. This variable is similar to 
the Bengtsson-Zerwas angle [T5], and is useful to separate g — > bb events because the 
radiated virtual gluon in the process Z° — > qqg is polarized in the plane of the three- 
parton event, and this is reflected in its subsequent splitting, by strongly favoring 
g — > qq emission out of this plane. Events with | cos 01234] > 0.8 are rejected. 

Thirdly, the b jets coming from a gluon tend to have lower energy than the other 
two jets in the event. Figure § shows the energy sum distribution of jets 1 and 2. We 
require the energy sum to be smaller than 36 GeV. 

Finally, c jets have lower Pj-corrected mass than b jets. Figure ^ shows the greater 
of the P^-corrected mass determined for jet 1 and jet 2 after the above cuts. Many 
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Source Efficiency (%) 

B 3.86 ±0.52 

C 0.10 ±0.02 

Q 0.73 ±0.05 

Table 1: Efficiencies after all cuts for the three categories. Errors are statistical only. 

background events lie below 2.0 GeV. Hence we require Mp T to be greater than 2.0 
GeV. 



6. Result 

After requiring all the above mentioned cuts, 62 events are selected in the data. The 
number of background events is estimated, using the Monte Carlo simulation, to be 
27.6, where 63% of the background comes from Z — > bb events, 27% from g — > cc 
events and the remaining 10% from Z —>■ qq (q ^ b) events. Table [| shows the tagging 
efficiencies for the three categories of events, where the errors are statistical only. From 
these efficiencies and the fraction of events selected in the data fd = (2.14±0.27) x 10~ 4 , 
the value of g hb can be determined: 

_ fd — (1 — 9cc)tQ - 9cctC /. \ 

9bb — • UJ 

e B - e Q 

We obtain 

g bi = (3.07 ± 0.71) xl0~ 3 , (2) 
where the error is statistical only. 



7. Systematic Errors 

The efficiencies for the three event categories are evaluated using the Monte-Carlo 
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simulation. The limitations of the simulation in estimating these efficiencies lead to an 
uncertainty on the result. The error due to the limited Monte-Carlo statistics in the 
efficiency evaluation is Ag bb = ±0.44 x 10~ 3 . This uncertainty comes mainly from the 
efficiency to tag Q events. 

A large fraction of events remaining after the selection cuts contain B and D 
hadrons. The uncertainty in the knowledge of the physical processes in the simula- 
tion of heavy-flavor production and decays constitutes a source of systematic error. 
All the physical simulation parameters are varied within their allowed experimental 
ranges. In particular, the h and c hadron lifetimes as well as production rates are var- 
ied, following the latest recommendations of the LEP Heavy Flavour Working Group 
|T5| |. The uncertainties are summarized in Table ^. 

The simulation of the signal events is based on the JETSET parton shower Monte 
Carlo, which is in good agreement with the theoretical predictions |IJ. In order to 
estimate the uncertainty on this assumption, we have produced 10,000 g — > bb events 



using GRC4F |L7| at the generator level. The signal tagging efficiency, e#, mainly 
depends on the energy of the gluon splitting into bb. This efficiency function, computed 
with JETSET, is reweighted by the ratio of the GRC4F to JETSET initial energy 
distributions to obtain the average efficiency. A systematic error of Ag bb = ±0.09 x 10~ 3 
is estimated from the difference in efficiency between the two Monte-Carlo models. 

The dependence of the efficiency on the 6-quark mass has also been investigated 
at the generator level. Events are generated using the GRC4F Monte Carlo, which 
is based on a matrix element calculation including 6-quark masses. The variation 
of €b was computed for 6-quark masses between 4.7 and 5.3 GeV/c 2 . The resultant 
uncertainty is estimated to be Ag bb = ±0.06 x 10~ 3 . 

The uncertainty in the ratio of the g cc background events, Ag c5 = ±0.40%, 
gives an error Ag bb = ±0.09 x 10 -3 . 

There is about a 5% discrepancy in the 4-jet rate between data and Monte Carlo 
at our y cut value. The uncertainty due to the discrepancy is estimated by increasing 
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the number of background events in the Monte Carlo and is found to be Ag bb = 
±0.14 x 1(T 3 . 

In the Monte-Carlo simulation charged tracks used in the topological vertex tag are 
smeared and rejected to reproduce better distributions in the data. Uncertainties in 
the efficiencies due to this smearing and rejection are assessed by evaluating the Monte- 
Carlo efficiencies without the smearing and rejection algorithm. The difference in the 
g bb result is taken as systematic error. The errors on g bb due to the tracking resolution 
and efficiency are then estimated to be Ag bb = ±0.26 x 10~ 3 and = ±0.29 x 10 -3 , 
respectively. 

Table |^ summarizes the different sources of systematic error on g bb . The total 
systematic error is estimated to be the sum in quadrature, 0.66 x 10~ 3 . 

8. Summary 

A preliminary measurement of the gluon splitting rate to a bb pair in hadronic 
Z° decays collected by SLD has been presented. Excellent SLC and VXD3 perfor- 
mance provides advantages not only for the 6-tag efficiency but also for the topological 
selections. The result is 

g bb = (3.07±0.71(stat.) ±0.66(syst.)) x 10" 3 (preliminary). 

where the first error is statistical and the second includes systematic effects. 
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Figure 1: Angular distribution between vertex axes in jet 1 and jet 2 (0.9 < cos #12) (points) 
The simulated distribution is shown clS ct histogram. 
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Figure 2: The distribution of the cosine of angle between the plane II12 formed by jets 1 
and 2 and the plane H34 formed by jets 3 and 4, for data (points) and Monte Carlo simulation 
(histogram). We reject | cos 01234] > 0.8 
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Figure 3: The distribution of the energy sum of jet 1 and jet 2. 
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Figure 4: Maximum P^-corrected mass distribution between jet 1 and jet 2 after jet-energy- 
sum cut. Points indicate data, open box signal, hatched boxes are backgrounds. 
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